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two-dimensional  external  disturbances  affects  the  mixing  field  in  plane  two-stream 
shear  layers.  Our  study  relates  to  the  effects  of  forcing  on  mixing  of  species,  the 
kind  of  mixing  that  is  necessary  for  chemical  reaction  and  combustion.  This  is  in 
contrast  to  many  previous  investigations  of  forced  shear  layers  which  have  concen¬ 
trated  on  the  effects  of  forcing  as  it  relates  to  the  momentum  transport  properties  of 
the  turbulent  region  [1-6].  These  properties  typically  include  the  layer  growth  rate, 
mean  and  rms  velocity  field,  and  turbulent  stresses. 

The  major  part  of  this  effort  involved  the  quantification  of  the  mixed-fluid  con¬ 
centration  field,  in  terms  of  its  probability  density  function  (pdf),  in  actively  forced 
shear  layers  using  laser  induced  fluorescence  (LIF)  diagnostics.  In  a  separate  smaller 
effort,  the  effect  of  curvature  on  the  2-D  structure  of  the  shear  layer  was  investi¬ 
gated.  While  the  main  goal  is  to  understand  the  basic  fluid  mechanical  phenomena 
in  each  case,  it  is  hoped  that  the  knowledge  gained  can  be  used  in  devising  active 
and  passive  means  to  modify/control  the  flow  and  the  mixing  fipld.  In  the  sections 
that  follow,  a  short  summary  of  our  main  results  to  date  are  presented.  The  details 
are  attached  as  Appendices  A  and  B. 


2.  SUMMARY  OF  RESULTS 
2.1  Composition  of  Mixed  Fluid  in  a  Forced  Shear  Layer 


To  document  the  effect  of  forcing  on  mixing,  a  series  of  experiments  were  per¬ 
formed  where  the  probability  density  function  (pdf)  of  the  concentration  field  was 
measured  under  different  forcing  conditions.  The  experiments  were  conducted  in  a 
gravity-driven  liquid  shear  layer  facility  where  forcing  was  achieved  by  oscillating 
one  of  the  two  free-stream  velocities  using  an  oscillating  bellows  mechanism  in  one 
of  the  supply  lines.  The  mixing  field  was  measured  using  the  laser  induced  fluores¬ 
cence  (LIF)  technique  in  the  passive  scalar  mode.  The  technique  is  essentially  tiie 
same  as  that  developed  previously  by  Koochesfahani  and  Dimotakis  [7],  but  extended 
to  two  dimensions.  Using  laser  sheet  illumination,  a  2-D  CCD  image  array,  and 
high  speed  digital  image  acquisition  techniques,  whole-field  concentration  data  were 
obtained  simultaneously  over  a  grid  of  512x480  points  in  a  plane  with  high  spatial 
and  temporal  resolution.  The  success  of  this  portion  of  our  work  was  in  part  due  to 
a  high-speed  digital  image  acquisition  system  obta’ricu  through  the  DURIP  equipment 
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grant  (AFOSR-89-0130). 

The  main  conclusions  from  a  shear  layer  forced  by  2-D  disturbances  are  sum¬ 
marized  below: 

•  Composition  distribution  of  mixed  fluid  in  forced  shear  layers  is  essentially  uni¬ 
form  across  the  width  of  the  shear  layer  similar  to  previous  results  in  natural 
layers. 

•  Forcing  changes  both  the  amount  and  the  composition  of  mixed  fluid.  For  the 
cases  studied,  forcing  results  in  an  increase  of  the  total  amount  of  mixed  fluid 
(integrated  across  the  layer  width)  and  a  shift  of  the  predominant  mixed-fluid 
concentration  to  larger  values. 

•  The  increase  in  the  amount  of  mixing  in  the  forced  layer  appears  to  be  mostly 
due  to  the  increase  in  the  layer  width  and  not  so  much  the  result  of  improved 
small-scale  mixing. 

•  Details  of  the  forcing  waveform  shape  have  a  significant  influence  on  the  struc¬ 
ture  of  the  flow  and  the  mixing  field.  Asymmetric  (i.e.  nonsinusoidal)  forcing 
can  result  in  a  much  larger  change  in  the  composition  field  than  that  due  to 
sinusoidal  forcing. 


The  details  of  these  findings  along  with  the  full  description  of  the  work  were 
presented  in  a  paper  at  the  lUTAM  Symposium  on  Fluid  Mechanics  of  Stirring  and 
Mixing  in  August  1990  in  La  Jolla,  California  [8].  This  paper,  attached  as  Appendix 
A,  has  been  accepted  for  publication  in  Physics  of  Fluids  A  and  will  appear  in  April 
1991.  This  work  also  forms  the  major  pan  of  the  M.S.  thesis  of  Mr.  Colin  MacKin¬ 
non. 

Our  study  of  shear  layer  forcing  with  asymmetric  waveform  has  recently  lead  to 
interesting  new  results.  We  have  reported  earlier  [8,  9]  that  asymmetric  forcing  with 
a  waveform  of  symmetry  5  =  70%  results  in  a  much  larger  growth  rate  and  mixed- 
fluid  composition  modification  than  sinusoidal  (S  =  50%)  forcing.  Asymmetric  forc¬ 
ing  with  5  =  30%,  however,  produces  a  "forced"  shear  layer  which  is  quite  similar 
to  a  "natural"  layer  [9].  In  an  ongoing  effort,  we  have  found  that  the  similarity 
applies  to  both  the  layer  growth  and  the  mixing  characteristics  (e.g.  amount  of  mixed 
fluid).  Upon  closer  examination,  it  was  noted  that  the  two  cases  of  5  =  30%  and 
70%  have  identical  frequency  contents  comprising  of  basically  a  fundamental  and  a 
subharmonic  (8  and  4  Hz  in  this  case).  The  main  difference  between  the  two 
waveforms  is  the  phase  between  the  fundamental  and  subharmonic;  the  fundamental- 
subharmonic  phase  for  5  =  30%  is  180  degrees  apart  from  that  for  5  =  70%.  The 
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importance  of  the  phase  between  fundamental  and  subharmonic  to  the  evolution  of 
the  shear  layer  beyond  the  splitter  plate  tip  has  been  known  for  some  time  [5,  10, 
11].  Most  of  these  studies  seem  to  have  concentrated  on  the  initial  instability  region 
of  the  shear  layer  where  the  fundamental  frequency  corresponds  to  the  initial  shear 
layer  most  amplified  frequency.  Our  results  support  the  notion  that  similar 
phenomena  may  occur  in  the  case  of  the  turbulent  shear  layer  farther  downstream.  In 
this  case,  the  fundamental  frequency  would  correspond  to  the  local  natural  (i.e.  vor¬ 
tex  passage)  frequency  of  the  turbulent  shear  layer.  Further  work  is  necessary  to 
clarify  the  details  of  the  processes  involved. 


2.2  Effect  of  Curvature  on  a  Two-Stream  Shear  Layer 

The  extensive  body  of  knowledge  which  exists  on  turbulent  flow  curvature 
effects  concentrates  on  the  characterization  of  the  velocity  field  (see,  for  example. 
Ref.  12)  whereas  much  less  is  known  about  the  mixing  field.  Besides  the  inherent 
interest  in  curvature  effects  on  fluid  mechanics  phenomena,  use  of  curvature  as  a 
possible  tool  for  mixing  control  needs  to  be  investigated.  Curvature  could  change 
the  mixing  field  by  modifying  the  overall  large  scale  (2-D)  entrainment  into  the  flow 
and  also  by  changes  in  the  3-D  structure  of  the  flow  due  to  the  Taylor-Goertler  insta¬ 
bility. 

We  performed  a  joint  theoretical  and  experimental  investigation  of  curvature 
effects  on  the  2-D  structure  of  the  shear  layer.  The  inviscid,  linear,  parallel-flow  sta¬ 
bility  analysis  of  spatially  growing  disturbances  was  utilized  to  study  the  instability 
characteristics  of  a  curved  shear  layer.  To  study  curvature  effects  experimentally, 
our  shear  layer  facility  was  modified  to  accommodate  a  curved  test  section.  Experi¬ 
ments  involved  LIF  flow  visuahzation  of  the  mixing  pattern  in  both  natural  and 
forced  curved  shear  layers.  Experiments  and  stability  analysis  lead  to  the  conclusion 
that  the  effect  of  curvature  (positive  or  negative)  on  the  2-D  structure  of  the  flow  is 
minimal  even  at  very  small  curvatures.  Details  of  this  work  were  reponed  in  an  ear¬ 
lier  progress  report  [9]  which  is  included  here  as  Appendix  B. 

Our  results  to  date  indicate  that  curvature  may  not  play  any  major  role  in  modi¬ 
fying  the  initial  development  of  the  shear  layer  vortices  and  the  fluid  entrained  into 
them.  The  influence  of  curvature  on  streamwise  vorticity  is  known  to  be  very 
strong,  however.  How  this  affects  the  molecular  mixing  field  requlics  a  nioie  exten¬ 
sive  investigation. 
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The  probability  density  function  (pdf)  of  mixed-fluid  composition  in  nonreacting 
forced  liquid  mixing  layers  is  measured  using  laser  induced  fluorescence  diagnostics. 
Results  show  that  the  composition  distribution  of  mixed  fluid  in  forced  shear  layers 
is  essentially  uniform  across  the  width  of  the  shear  layer  similar  to  previous  results 
in  natural  shear  layers.  For  the  cases  studied,  forcing  results  in  an  increase  of  the 
total  amount  of  mixed  fluid  (integrated  across  the  layer  width)  and  a  shift  of  the 
predominant  mixed-fluid  concentration  to  larger  values.  The  increase  of  the  amount 
of  mixing  in  the  forced  layer  appears  to  be  mostly  due  to  the  increase  in  the  layer 
width  and  not  so  much  the  result  of  improved  small-scale  mixing.  It  is  shown  that 
the  details  of  the  forcing  waveform  shape  can  have  a  significant  influence  on  the 
structure  of  the  flow  and  the  mixing  field. 
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I.  INTRODUCTION 


It  is  well  documented  that  forcing  can  be  used  to  manipulate  the  evolution  of 
turbulent  mixing  layers.  Using  a  variety  of  techniques,  it  has  been  shown,  for  exam¬ 
ple,  that  two-dimensional  perturbations  may  suppress  or  enhance  the  spreading  rate  of 
the  shear  layer  [1-6].  Except  for  a  few  cases,  which  will  be  discussed  shortly,  virtu¬ 
ally  all  experimental  studies  of  forced  shear  layers  to  date  have  concentrated  on  the 
effects  of  forcing  on  mixing  as  it  relates  to  the  momentum  transpon  properties  of  the 
turbulent  region.  These  properties  typically  include  the  layer  growth  rate,  mean  and 
rms  velocity  field,  and  turbulent  stresses.  The  study  presented  here  relates  to  the 
effects  of  forcing  on  mixing  of  species,  the  kind  of  mixing  that  is  necessary  for 
chemical  reactions  and  combustion. 

We  are  aware  of  only  two  experimental  studies  where  information  is  available 
on  molecular  mixing  of  species  in  a  two-dimensionally  forced  shear  layer.  In  a  work 
by  Wygnanski,  et  al.  [1],  heating  was  used  as  a  passive  contaminant  mainly  for  the 
purpose  of  observing  the  structure  of  the  flow.  Average  and  rms  temperature  profiles 
were  used  to  conclude  that,  in  their  experiment,  the  forced  shear  layer  consisted  of 
an  orderly  array  of  vortices  whose  core  was  fairly  well  mixed.  The  work  of  Roberts 
[7,  8]  is  the  most  extensive  body  of  research  on  mixing  in  forced  shear  layers. 
Using  chemically  reacting  techniques  in  a  liquid-phase  shear  layer,  he  showed  that 
the  extent  of  molecular  mixing,  as  manifested  by  the  total  amount  of  chemical  pro¬ 
duct  in  the  layer,  can  be  altered  upon  shear  layer  forcing.  All  of  these  measurements 
were  performed  at  one  fixed  value  of  the  reactant  stoichiometric  ratio.  It  should  be 
noted,  however,  that  it  is  nearly  impossible  to  obtain  unambigous  information  on 
molecular  mixing  from  a  chemically  reacting  experiment  at  a  single  stoichiometric 
ratio.  The  source  of  the  difficulty  is  that  the  amount  of  chemical  product  depends  on 
both  the  amount  and  composition  of  the  molecularly-mixed  fluid.  A  full  description 
of  the  effect  of  forcing  on  molecular  mixing  requires  information  on  the  entire  com¬ 
position  distribution  of  the  mixed  fluid,  namely  the  probability  density  function  (pdf) 
of  the  concentration  field.  It  is  the  purpose  of  this  study  to  provide  such  informa¬ 
tion.  We  draw  attention  to  the  utility  of  this  type  of  information  in  mixing  control 
where  manipulating  a  turbulent  shear  layer  to  have  more  mixed  fluid  may  prove 
ineffective  if  the  resulting  mixed  fluid  has  an  undesirable  composition. 

In  the  current  work,  the  pdf  of  the  concentration  field  in  a  forced  liquid  shear 
layer  is  estimated  using  a  passive  scalar  technique.  It  is  known  that  this  technique 
yields  an  upper  bound  to  the  actual  molecular  mixing  (e.g.  see  Refs.  10,  11).  In 
using  this  technique,  a  passive  scalar  contaminant,  such  as  dye,  is  premixed  with  one 
of  the  streams.  The  concentration  of  dye  within  a  small  sampling  volume. 
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determined  by  the  spatial  and  temporal  resolution  characteristics  of  the  measuring 
apparatus,  is  then  recorded  as  a  function  of  time.  The  difficulty  arises  if  the  sam¬ 
pling  volume  is  larger  than  the  smallest  mixing  (diffusion)  scale,  as  is  usually  the 
case  in  high  Reynolds  number  flows.  Under  these  circumstances,  it  is  impossible  to 
determine  whether  the  two  fluids  are  mixed  or  not  within  the  measurement  resolu¬ 
tion.  For  the  measurements  reported  here,  for  example,  the  spatial  resolution  was 
estimated  to  be  approximately  4.4  times  larger  than  the  smallest  diffusion  scale  (see 
Section  2).  The  problem  of  finite  sampling  volume  can,  of  course,  be  solved  by 
using  diagnostics  that  rely  on  chemically  reacting  techniques  [10,  11]. 

With  all  the  potential  drawbacks  associated  with  the  passive  scalar  method  in 
determining  the  true  extent  of  molecular  mixing,  it  nevertheless  is  the  simplest  way 
to  get  an  estimate  of  the  pdf  of  the  composition  field  in  a  single  experiment.  While 
it  is  possible  to  determine  the  pdf  from  chemical-reaction  data  with  much  better 
accuracy  in  representing  molecular  mixing,  this  would  require  many  experiments  at 
different  reaction  stoichiometric  ratios.  The  results  presented  in  this  paper  are  not 
intended  to  provide  an  absolute  measure  of  the  extent  of  mixing;  as  mentioned  ear¬ 
lier,  they  yield  an  upper  bound  to  the  actual  molecular  mixing.  The  primary  focus 
of  these  results  is  to  characterize  the  relative  changes  in  the  pdf  of  the  composition 
field  in  a  forced  shear  layer  compared  to  a  base  state  (the  natural  layer,  in  this  case). 
We  believe  such  comparisons  are  warranted  since  the  relative  resolution  (compared  to 
the  diffusion  scale)  is  essentially  the  same  in  all  of  these  measurements. 

The  paper  is  organized  as  follows.  Section  2  gives  a  brief  description  of  the 
measuring  technique  and  instrumentation.  Discussion  of  results  in  Section  3  includes 
both  sinusoidal  and  non-sinusoidal  forms  of  forcing.  Most  of  the  quantities  discussed 
are  computed  from  the  measured  pdf’s.  The  definitions  of  these  quantities  are  given 
in  the  Appendix. 
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n.  EXPERIMENTAL  FACILITY  &  INSTRUMENTATION 


Experiments  were  conducted  in  a  gravity-driven  water  shear  layer  apparatus  in 
which  the  free-stream  fluids  were  supplied  from  two  independent  reservoirs  (see  Fig¬ 
ure  1).  The  test  section  had  a  cross  section  of  4  cm  (height)  x  8  cm  (span)  and  was 
35  cm  long.  Forcing  was  achieved  by  oscillating  the  high-speed  free-stream  speed 
using  an  oscillating  bellows  in  the  high-speed  stream  supply  line  as  shown  in  Figure 
2.  The  motion  of  the  bellows  was  controlled  by  a  shaker  coil  mechanism  whose 
command  signal  originated  from  a  function  generator  (HP3314A). 

For  the  results  described  here,  the  free-stream  speeds  were  set  to  t/j  =  40  cm/s 
and  (72  =  20  cm/s.  The  natural  shear  layer  roll-up  frequency  F^,  under  these  condi¬ 
tions  was  about  27  Hz.  The  bellows  was  driven  sinusoidally  at  frequencies  of 
F  =  4,  8  Hz.  In  the  case  of  4  Hz,  a  skewed  sine  wave  motion  was  also  used  in 
order  to  observe  the  effect  of  nonsinusoidal  forcing.  In  this  case,  we  characterize  the 
bellows  motion  in  terms  of  the  symmetry  parameter  S,  which  is  the  percentage  of  a 
period  (in  one  cycle)  during  which  the  bellows  is  moving  in  the  downward  (i.e. 
increasing  free-stream  speed)  direction.  Results  are  presented  for  S  =  70%  (F  =  4 
Hz).  Note  that  S  =  50%  corresponds  to  sinusoidal  forcing.  The  rms  velocity 
fluctuation  imposed  on  the  high-speed  stream  as  a  result  of  forcing  was  estimated  to 
be  <  0.5%. 

The  mixing  field  was  measured  using  the  laser  induced  fluorescence  (LIF)  tech¬ 
nique  described  by  Koochesfahani  &  Dimotakis  [9,  11].  A  fluorescent  dye  (disodium 
fluorescein)  was  premixed  with  the  low-speed  free-stream  fluid  which  was  subse¬ 
quently  diluted  in  the  shear  layer  as  a  result  of  mixing  with  the  fluid  from  the  high¬ 
speed  free  stream.  Recording  the  fluorescence  intensity  allowed  a  quantitative  meas¬ 
urement  of  the  dye  concentration  and,  therefore,  the  relative  concentration  of  high¬ 
speed  to  low-speed  fluid  in  the  layer.  This  was  done  by  noting  that  the  local  instan¬ 
taneous  dye  concentration  in  a  sampling  volume  is  given  by 


^2 


"Ui  -I-  1)2 


where  is  the  free-stream  dye  concentration  carried  by  the  fluid  on  the  low-speed 
side  and  \)i  and  U2  are  the  volumes  of  fluid  from  the  high-speed  and  low-speed 
streams,  respectively,  in  the  sampling  volume.  The  normalized  concentration  of 
iiigh-speed  fluid  ^  is  thus  equal  to 


-4- 


4=  1-Q/c,, 


where  ^  =  'Ui/('Ui +1)2)  is  defined  as  the  high-speed  fluid  volume  fraction. 

The  LIF  measurements  were  carried  out  over  a  plane  defined  by  a  laser  sheet 
aligned  along  the  flow  direction  at  the  shear  layer  midspan  location.  The  laser  sheet 
(thickness  about  0.5  mm)  was  formed  by  passing  the  beam  of  a  4  W  argon-ion  laser 
(EXCEL-3000)  through  a  combination  of  spherical  and  cylindrical  lenses.  The 
fluorescence  intensity  was  recorded  by  a  2-D  CCD  camera  (NEC  TI-24A)  operating 
at  60  fields/s  with  an  exposure  of  2  msec.  The  4  cm  width  of  the  test  section  was 
imaged  on  200  pixels  (out  of  480  vertical  pixels  available)  of  the  camera  so  that  the 
spatial  resolution  was  200  |im  x  200  |J.m.  At  this  resolution,  the  512  horizontal  pix¬ 
els  of  the  camera  allowed  measurements  along  10.24  cm  of  the  flow  in  the  stream- 
wise  direction. 

The  camera  output  was  digitized  to  8  bits  at  full  video  rate  into  hard  disk  in 
real  time  by  a  digital  image  acquisition  system  (TRAPDC-5500).  Note  that,  with  this 
arrangement,  simultaneous  multipoint  measurements  of  species  concentration  were 
obtained  at  512  x  480  grid  points  (512  x  200  pixels  inside  test  section)  as  a  function 
of  time  with  high  spatial  (200  |ini  x  200  temporal  (2  msec)  resolution.  For 

each  run,  256  sequential  LIF  images  (i.e.  64  MBytes  of  data)  were  acquired  to  con¬ 
struct  the  pdf  and  also  allow  the  monitoring  of  the  spatiotemporal  evolution  of  the 
concentration  field  in  two  dimensions.  The  number  of  structures  passing  the  field  of 
view  during  the  data  acquisition  period  was  approximately  34  for  the  case  of  forcing 
at  4  Hz  and  68  for  forcing  at  8  Hz. 

The  steps  for  processing  the  fluorescence  intensity  data  are  described  in  detail 
by  Koochesfahani  &  Dimotakis  [9]  and  are  not  repeated  here.  We  should  mention 
that  in  the  present  experiments  the  laser  intensity  attenuation  due  to  absorption  by  the 
dye  was  negligible;  the  non-uniformity  in  the  laser  sheet  illumination  and  camera 
pixel  sensitivity  were  taken  into  account,  however. 

For  the  measurements  presented  here,  the  region  of  the  flow  between 
15  cm  <  X  <  25  cm  was  imaged  on  the  camera.  At  the  middle  of  this  region  (i.e. 
at  X  =20  cm ),  the  local  natural  shear  layer  width  5i  (see  Appendix)  was  estimated 
to  be  about  3.3  cm,  resulting  in  a  local  Reynolds  number  of  6,600.  This  Reynolds 
number  corresponds  to  a  natural  shear  layer  during  the  mixing  transition  [10]. 
Recent  results  by  Miller  and  Dimotakis  [12]  indicate  that  the  smallest  expected  scalar 
diffusion  scale  (Batchelor  scale)  Xq  is  very  close  to  the  Kolmogorov  scale  X^  in 
water  where  the  Schmidt  number  is  approximately  600.  Schmidt  number.  Sc  =  v/D , 
is  the  ratio  of  the  diffusion  coefficient  of  momentum  to  that  of  mass.  Estimating  the 
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_3^ 

Kolmogorov  scale  by  =  Sj  Re  ^ ,  gives  a  value  for  the  smallest  diffusion  scale 
Xq  ~  45  p.m  for  the  conditions  at  a:  =20  cm  mentioned  above.  As  discussed  in  the 
introduction,  our  results  provide  an  upper  bound  to  the  actual  molecular  mixing  since 
the  typical  spatial  resolution  of  these  measurements  is  about  4.4  times  larger  than  the 
smallest  mixing  scale.  The  characterization  of  the  relative  changes  in  the  pdf  of  the 
composition  field  in  a  forced  shear  layer  compared  to  the  natural  layer  is  expected  to 
be  reliable,  however,  since  the  relative  resolution  (compared  to  the  diffusion  scale)  is 
basically  the  same  in  all  of  these  measurements. 
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ni.  RESULTS  &  DISCUSSION 

Typical  digital  planar  LIF  data  for  natural  and  forced  shear  layers  are  presented 
in  Figures  3  through  6.  Each  figure  shows  a  time  sequence  (Ar  =  1/15  sec)  of  the 
mixing  field  over  a  downstream  region  15  cm  <  x  <  25  cm .  These  time  sequences 
represent  less  than  1%  of  the  data  available  for  a  given  run.  Flow  is  from  right  to 
left  with  the  high-speed  stream  on  top.  At  each  time  step,  the  image  is  actually 
composed  of  50,000  (200  in  y  and  500  in  x )  pixels.  Colors  are  assigned  to  different 
mixed  fluid  concentration  levels  in  order  to  aid  the  interpretation  of  data.  The  color 
assignments  are  as  follows: 


0.88  <  ^  <  0.94 

Red 

0.64  <  ^  <  0.70 

Dark  green 

0.82  <  ^  <  0.88 

Orange 

0.58  <  ^  <  0.64 

Purple 

0.76  <  ^  <  0.82 

Yellow 

0.52  <  ^  <  0.58 

Blue 

0.70  <  ^  <  0.76 

Green 

0.40  <  ^  <  0.52 

White 

0.28  <  ^  <  0.40 

Crimson 

Concentration  levels  below  0.28  were  not  labeled  since  the  probability  of  their 
occurence  was  relatively  low.  Note  also  that  pure  fluids  from  the  high-speed  (q  =  1) 
and  low-speed  (q  =  0)  streams  were  not  assigned  colors  and  appear  black. 

Figures  3-5  show  the  enhanced  growth  rate  of  the  sinusoidally  forced  shear 
layer  compared  to  the  natural  case.  This  observation  is  consistent  with  previous 
results  [3,4,6]  establishing  the  large  increase  of  layer  growth  rate  under  low- 
frequency  (relative  to  natural  frequency  F^)  forcing.  Additionally,  present  results 
indicate  that  using  a  nonsinusoidal  form  of  forcing  can  lead  to  even  larger  growth 
rates;  compare  Figures  5  and  6.  In  the  case  of  sinusoidal  forcing  (S  =  50%)  at  4 
and  8  Hz,  we  see  a  single  large  scale  vortex  lump  passing  through  the  middle  of  the 
field  of  view  (;c  =  20  cm).  For  the  particular  form  of  nonsinusoidal  forcing  used 

here  (F  =  4  Hz,  S  =  70%),  the  flow  is  characterized  by  large  pairing  vortices  con- 

vecting  through  the  field  of  view  (Fig.  6).  Note  that  the  size  of  the  structures  in 
these  forced  layers  has  become  comparable  to  the  test  section  height.  It  is  quite  pos¬ 
sible  that  the  finite  height  of  the  test  section  may  be  interfering  with  the  growth  and 

evolution  of  the  structures.  This  is  particularly  serious  in  the  case  of 
F  =  4  Hz,  S  =  70%  where  the  bottom  side  of  the  structures  actually  comes  in  con¬ 
tact  with  the  lower  wall  of  the  test  section  at  some  downstream  distance. 

From  the  image  time  series  discussed  above,  the  probability  density  function 
(pdf)  of  the  concentration  field  can  be  constructed  at  any  of  the  50,000  measurement 
points  in  the  (x,  y)  plane.  For  the  remainder  of  this  paper,  we  will  discuss  the 
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results  only  at  x  =20  cm  (i.e.  midway  between  the  left  and  right  edges  of  the  LIF 
images  in  Figures  3-6).  These  results  can  be  considered  representative  because  of 
the  small  variation  of  the  pdf’s  over  the  range  of  streamwise  locations  investigated 
here.  Note  that  the  descriptions  of  the  quantities  discussed  below  and  how  they  are 
computed  from  the  pdf’s  are  included  in  the  Appendix. 

The  composition  distribution  for  natural  and  forced  (F  =  4  Hz)  shear  layers  are 
shown  in  Figure  7(a-c)  in  terms  of  the  contour  plot  of  the  pdf,  p(^,  y),  at  x  =20 
cm.  One  main  feature  is  apparent  in  this  figure:  in  all  cases,  the  composition  distri¬ 
bution  of  mixed  fluid  is  uniform  across  the  width  of  the  shear  layer.  This  is  also 
true  for  the  case  of  F  =  8  Hz  whose  pdf  was  not  included  in  Figure  7.  The  unifor¬ 
mity  of  mixed-fluid  composition  in  natural  shear  layers,  known  for  some  time  now 
(e.g.  see  Refs.  9,  11,  16),  is  an  attribute  of  the  vortical  motion  of  the  large  scale 
structures.  That  similar  results  also  apply  to  forced  shear  layers  is,  therefore,  not 
surprising.  These  findings  are  consistent  with  the  recent  direct  numerical  simulation 
of  mixing  in  a  forced  shear  layer  by  Buell  &  Mansour  [13]. 

One  advantage  of  the  uniformity  of  mixed-fluid  pdf  is  that  the  composition  field 
of  the  whole  layer  can  be  characterized  by  a  total  pdf,  F(^),  which  is  the  integral  of 
p(^,y)  across  the  width  of  the  layer  (see  Appendix).  The  total  pdfs,  presented  in 
Figure  8,  illustrate  the  main  effects  of  forcing  on  the  mixing  field.  Note  that  forcing 
changes  both  the  amount  and  the  composition  of  the  mixed  fluid.  For  the  cases  we 
have  studied  thus  far,  the  predominant  concentration  shifts  to  larger  values,  namely 
mixing  is  enhanced  preferentially  at  higher  values  of  Visual  inspection  of  the 
total  pdf’s  suggests  that  the  total  amount  of  mixed  fluid  in  the  layer  has  increased 
upon  forcing;  quantitative  measure  will  be  presented  later.  An  interesting  result  is 
that  the  details  of  the  forcing  waveform  shape  have  a  significant  influence  on  the 
structure  of  the  flow,  as  discussed  earlier,  and  the  mixing  field.  Figure  8b  shows 
that  asymmetric  forcing  (i.e.  nonsinusoidal)  can  result  in  a  much  larger  change  in  the 
composition  field  than  that  due  to  sinusoidal  forcing.  This  can  also  be  visually 
confirmed  by  observing  the  increased  amount  of  high  concentration  fluid  (i.e.  red, 
orange,  yellow  colors)  in  Figure  6  compared  to  Figure  3. 

As  was  discussed  above,  in  all  forcing  cases  presented  here  the  predominant 
mixed-fluid  composition  shifts  to  larger  values.  This  may  be  interpreted  to  be  due  to 
an  increased  asymmetry  of  the  entrainment  ratio  into  the  shear  layer.  It  is  known 
that  the  natural  shear  layer  entrains  an  unequal  amount  of  fluid  from  the  two  free 
streams.  The  entrainment  ratio  E  is  defined  to  be  the  volume  ratio  of  the  entrained 
high-speed  to  low-speed  fluids.  In  the  recent  model  of  entrainment  by  Dimotakis 
[14],  the  geometric  properties  of  the  large-scale  flow  structures  of  the  natural  mixing 
layer  were  utilized  to  derive  an  expression  for  the  entrainment  ratio.  For  the  case  of 
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uniform  density,  the  expression  was  given  as  £  =  1  +  3.9  djx,  where  is  the  vor- 
ticity  thickness.  This  expression  connects  the  entrainment  ratio  directly  to  the  layer 
growth  rate.  The  equivalent  relationship  applicable  to  forced  layers  has  not  been 
established  to  our  knowledge  and  whether  the  expression  for  natural  shear  layer 
entrainment  ratio  can  be  applied  to  forced  layers  remains  to  be  proven.  We  note, 
however,  that  compared  to  the  natural  layer,  the  forced  shear  layers  studied  here  have 
larger  growth  rates;  they  also  have  higher  values  of  predominant  mixed-fluid  concen¬ 
tration. 

The  pdf  of  the  composition  field  p(^,  y)  can  be  used  to  calculate  the  transverse 
profiles  of  some  quantities  of  interest  (see  Appendix).  The  reference  location  y  =  0 
and  local  layer  thickness  8i  used  in  these  profiles  were  determined  according  to  the 
description  in  the  Appendix.  The  average  concentration  ^(y)  is  shown  in  Figure  9. 
The  profiles  for  forced  conditions  develop  a  "flat"  region  similar  to  that  observed  by 
Wygnanski,  et  al.  [1],  who  attributed  this  feature  to  an  orderly  array  of  vortices  with 
fairly  well  mixed  cores.  The  probabilities  of  finding  unmixed  fluid  from  low-  and 
high-speed  streams,  p^iy)  and  Pi(y)  respectively,  are  displayed  in  Figure  10.  Note 
the  large  increase  in  the  probability  of  finding  unmixed  fluid  (from  both  streams)  in 
the  middle  of  the  layer  as  forcing  is  applied.  Examining  the  LIF  images  in  Figures 
3  through  6  leads  to  the  conclusion  that  the  increase  is  due  to  the  more  pronounced 
and  larger  entrainment  tongues  between  the  large  structures  in  the  forced  layer. 
Correspondingly,  the  total  probability  of  finding  mixed  fluid  p^iy)  in  Figure  11 
shows  a  decrease  in  the  middle  of  the  forced  shear  layer,  indicating  a  reduction  in 
the  amount  of  mixed  fluid. 

As  a  measure  of  the  total  amount  of  mixed  fluid  across  the  layer,  the  mixed- 
fluid  thickness  5^  was  calculated.  This  measure  is  essentially  the  area  under  the 
Pm(y)  curve.  The  table  below  compares  the  total  amount  of  mixed  fluid  across  the 
layer  and  also  the  amount  per  unit  width  of  the  layer  for  the  cases  investigated. 


6„(mm) 

5./ 

Natural 

16.1 

0.48 

F  =  8  Hz,  5  =  50% 

19.3 

0.51 

F  =  4  Hz,  5  =  50% 

18.3 

0.45 

F  =  4  Hz,  5  =  70% 

19.7 

0.43 

The  most  imponant  result  is  that  while  the  total  amount  of  mixed  fluid  in  the  forced 
layer  has  increased  by  as  much  as  22%,  the  amount  of  mixing  per  unit  width  of  the 
layer  has  remained  nearly  constant.  This  suggests  that  the  forced  layers  in  our  study 
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have  more  mixed  fluid  mostly  by  virtue  of  the  increased  local  width  of  the  layer  and 
not  so  much  because  of  improved  small  scale  mixing.  In  fact,  as  noted  earlier,  there 
is  less  mixed  fluid  in  the  center  of  a  forced  shear  layer  compared  to  the  natural 
layer. 

The  effect  of  2-D  forcing  on  the  development  of  3-D  small  scales  need  to  be 
studied  in  greater  detail.  Recent  results  by  Huang  &  Ho  [15]  connect  the  production 
of  random  small-scale  eddies  to  the  interaction  between  the  merging  spanwise  struc¬ 
tures  and  the  streamwise  vortices.  The  present  results  suggest  that  2-D  forcing  may 
not  lead  to  any  increase  of  small  scale  mixing.  It  is  not  clear  whether  this  is  true 
for  all  forced  shear  layers.  In  particular,  we  do  not  have  any  information  on  the 
effects  of  forcing  amplitude  on  the  mixing  field.  Furthermore,  much  work  still 
remains  to  understand  how  the  frequency  contents  (i.e.  shape)  of  forcing  waveform 
affect  the  large  and  small  scales  in  a  turbulent  shear  layer. 
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rv.  CONCLUSIONS 


The  results  presented  show  that  the  composition  distribution  of  mixed  fluid  in 
forced  shear  layers  is  essentially  uniform  across  the  width  of  the  shear  layer  similar 
to  previous  results  in  natural  shear  layers.  For  the  cases  studied,  the  total  amount  of 
mixed  fluid  integrated  across  the  layer  width  increases.  The  amount  of  mixed  fluid 
per  unit  width  of  the  shear  layer  does  not  vary  significantly,  however.  This  implies 
that  the  increase  of  mixing  in  the  forced  layer  is  mostly  due  to  the  increase  in  the 
layer  width  and  not  because  of  improved  small-scale  mixing.  In  all  the  cases  inves¬ 
tigated  here,  mixed-fluid  composition  in  the  shear  layer  changes  with  forcing  such 
that  the  predontinant  mixed-fluid  concentration  shifts  to  larger  values.  In  other 
words,  mixing  is  enhanced  preferentially  at  higher  values  of  concentration.  Finally, 
results  show  that  nonsinusoidal  forcing  can  have  a  significant  influence  on  the  struc¬ 
ture  of  the  flow  and  the  mixing  field. 
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APPENDIX 


This  section  gives  the  definitions  of  the  quantities  discussed  in  this  paper  and 
how  they  are  related  to  the  pdf  of  the  high-speed  fluid  concentration  Most  of 
these  definitions  are  based  on  those  described  by  Koochesfahani  &  Dimotakis  [9,  11]. 

The  pdf  of  ^  at  a  given  point  y  in  the  shear  layer  is  denoted  by  p(^,y).  Since 
the  pdf  is  normalized,  we  have 


1 

J  pi^,y)d^=\.  (Al) 

0 

The  average  concentration  of  high-speed  fluid  ^  is  then  given  by 

?Cy)  =  J  (A2) 

0 

In  the  ideal  case,  the  two  ’delta’  functions  in  p{^,y)  at  ^  =  0,  1  are  associated 
with  the  pure  unmixed  fluids  from  the  low-  and  high-speed  sides,  respectively.  In 
practice,  however,  the  width  e  of  these  delta  functions  is  nonzero  and  is  dictated  by 
the  overall  signal-to-noise  ratio  of  the  measurement.  In  the  present  work,  the  range 
of  ^  was  divided  into  21  levels  so  that  the  value  of  e  was  about  0.047.  Concentra¬ 
tions  in  the  range  0<^<e  are  assigned  to  pure  (unmixed)  low-speed  fluid  and 
l-e<^<l  to  pure  high-speed  fluid.  The  range  E<^<l-e  would,  therefore 
correspond  to  mixed  fluid.  With  this  description,  the  probabilities  p^  and  p^  of 
finding  unmixed  fluids  from  low-  and  high-speed  sides  respectively,  are  given  by 

£  1 

Po(y)  =  J  p(^>y)d^ ,  Pi(y)  =  f  p(^,y)dt  (A3) 

0  1  -£ 

The  mixed-fluid  pdf  corresponds  to  the  ponion  of  p(^,y)  excluding  the  two 
delta  functions.  The  area  under  the  mixed-fluid  pdf,  p„,  gives  the  total  probability 
of  finding  mixed  fluid  at  any  concentration.  It  is  given  by 

1  -e 

Pm(y)  =  j  Pi%^y)d%.  (A4) 

€ 

A  value  of  p^iy)  of  less  than  unity  implies  the  presence  of  unmixed  fluid.  The 
reference  location  y  =  0  which  is  used  throughout  this  work  was  arbitrarily  chosen  to 
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coincide  with  the  point  at  maximum  Pm(y)-  Additionally,  the  reference  thickness  Sj 
which  is  essentially  equivalent  to  the  layer  visual  thickness  (see  Ref.  11)  was  defined 
as  the  1%  width  of  p^iy)  where  p„  has  dropped  to  1%  of  its  maximum  value. 

Since  it  was  observed  that  the  shape  of  the  mixed-fluid  pdf  was  nearly  invariant 
across  the  width  of  the  shear  layer  for  all  the  cases  studied,  the  composition  field  of 
the  whole  layer  can  be  characterized  by  a  single  pdf.  The  total  pdf  integrated  across 
the  layer  is  given  by 

=  J  pi^,y)dy,  (A5) 

The  total  amount  of  mixed-fluid  integrated  across  the  shear  layer  width  can  be 
characterized  by  the  mixed-fluid  thickness  5^  defined  by 

5m  =  J  Pm(y)dy=j^  (A6) 

In  the  case  of  chemically  reacting  shear  layers,  an  often-used  measure  of  the  total 
amount  of  chemical  product  in  the  layer  is  the  product  thickness.  For  fast  chemical 
reactions  in  the  limit  of  large  and  small  reaction  stoichiometric  ratios,  two  product 
thicknesses  are  defined  as  follows  (see  Refs.  11,  16) 

=  f  I  %P(^y)dkdy,  5p^  =  |  J  {\-^)p{\,y)d'^  dy.  (A7) 

-oo  £  -oo  £ 

It  is  easy  to  show  that  the  mixed-fluid  thickness  and  the  two  product  thicknesses  are 
related  by 

8".  =  +  ht  <A8) 

We  note  that  it  is  the  product  thickness  5^^  which  Robens  [7,  8]  uses  as  the  variable 
to  describe  the  effect  of  forcing  in  his  forced  shear  layer  studies. 
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FIGURE  CAPTIONS 


Figure  1.  Schematic  of  the  mixing  layer  facility. 

Figure  2.  Shear  layer  forcing  mechanism. 

Figure  3.  Digital  LIF  image  time  sequence  of  the  mixing  field  for  the  natural  layer. 
Flow  is  from  right  to  left  and  time  increases  from  top  to  bottom  (At  =  1/15  sec). 

Figure  4.  Digital  LIF  image  time  sequence  of  the  mixing  field  for  sinusoidal  forcing 
at  F  =  8  Hz,  S  =  50%. 

Figure  5.  Digital  LIF  image  time  sequence  of  the  mixing  field  for  sinusoidal  forcing 
at  F  =  4  Hz,  5  =  50%. 

Figure  6.  Digital  LIF  image  time  sequence  of  the  mixing  field  for  nonsinusoidal 
forcing  at  F  =  4  Hz,  S  =  70%. 

Figure  7.  Contour  plot  of  the  composition  distribution  p(^,  y)  at  x  =20  cm, 

(a)  Natural;  (b)  Forced,  F  =  4  Hz,  sinusoidal  (S  =  50%); 

(c)  Forced,  F  =  4  Hz,  nonsinusoidal  (S  =  70%). 

Figure  8.  Effect  of  forcing  on  the  total  pdf  of  the  concentration  field, 

(a)  Sinusoidal  forcing;  (b)  Nonsinusoidal  forcing. 

Figure  9.  Transverse  profile  of  average  high-speed  fluid  concentration. 

Figure  10.  Transverse  profiles  of  pure  free-stream  fluid  probability, 

(a)  Low-speed  fluid;  (b)  High-speed  fluid. 

Figure  11.  Transverse  profile  of  total  mixed-fluid  probability. 
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1.  INFLUENCE  OF  FORCING  ON  THE  COMPOSITION  OF 
MIXED  FLUID  IN  A  TWO-STREAM  SHEAR  LAYER 


It  is  well  documented  that  forcing  can  greatly  modify  the  behavior  of  turbulent 
plane  mixing  layers  [1-6].  While  most  of  the  research  in  this  area  has  concentrated 
on  the  effects  of  forcing  on  the  velocity  field  and  layer  growth  rate,  the  modification 
of  the  molecular  mixing  field  has  received  much  less  attention.  Roberts  [4,5]  has 
shown  that  the  extent  of  molecular  mixing,  as  manifested  by  the  amount  of  chemical 
product  formed  in  the  layer,  can  be  altered  by  shear  layer  excitation.  It  is  not 
known,  at  this  time,  how  the  composition  distribution  of  the  molecularly-mixed  fluid 
in  a  shear  layer  is  affected  upon  forcing.  It  is  important  to  note  that  manipulating  a 
turbulent  shear  layer  to  have  more  mixed  fluid  may  be  ineffective  if  the  resulting 
mixed  fluid  has  an  undesirable  composition.  It  is  the  purpose  of  this  study  to  deter¬ 
mine  how  the  probability  density  function  ipdf)  of  the  concentration  field  is  affected 
when  a  turbulent  shear  layer  is  forced.  The  investigation  naturally  addresses  the 
behavior  of  both  the  large  scale  (2-D)  entrainment  and  the  small  scale  (3-D)  structure 
of  the  forced  flow  since  they  determine  the  distribution  of  mixed  fluid  in  the  layer. 

We  are  conducting  experiments  in  a  gravity-driven  liquid  shear  layer  facility 
where  forcing  is  achieved  by  oscillating  one  of  the  two  free-stream  velocities  using 
an  oscillating  bellows  mechanism  in  one  of  the  supply  lines  [7].  Utilizing  laser 
induced  fluorescence  (LIF)  diagnostics  [8],  we  monitor  the  relative  concentration  of 
high-speed  to  low-speed  fluid  in  the  shear  layer  in  terms  of  high-speed  fluid  volume 
fraction  For  the  results  described  here,  the  shear  free-stream  speeds  were  set  to 
Uj  =  40  cm/s,  U2  =  20  cm/s.  The  fluorescence  intensity  over  a  plane  illuminated  by 
a  laser  sheet  was  recorded  by  a  2-D,  CCD  array  with  a  spatial  resolution  of  200  |im 
and  temporal  resolution  of  2  msec  at  60  fields/s.  The  array  output  was  digitized  to  8 
bits  at  full  video  rate  into  hard  disk  in  real  time  by  a  digital  image  acquisition  sys¬ 
tem  (TRAPIX-5500).  Typically  240  sequential  LIF  images  (i.e.  60  MBytes  of  data) 
were  acquired  to  ensure  the  reliability  of  statistics  and  also  allow  the  monitoring  of 
the  spatiotemporal  evolution  of  the  concentration  field  in  two  dimensions.  These 
digital  images  provide  spatially  and  temporally  resolved  concentration  data  at  about 
50,000  (simultaneous)  points  over  a  plane  in  the  shear  layer  as  a  function  of  time. 

We  present  in  Figures  1  through  3  the  effect  of  forcing  on  the  shear  layer  struc¬ 
ture.  Each  figure  shows  a  time  sequence  (Ar  =  1/15  sec)  of  the  mixing  field  over  a 
downstream  region  16.7  cm  <  x  <  24.5  cm.  At  each  time  step,  the  image  is  actu¬ 
ally  composed  of  40,000  (100  in  y  and  400  in  x)  point  measurements  of  concentra¬ 
tion  which  have  been  used  to  calculate  a  simulated  fast  chemical  reaction  at  unity 
equivalence  ratio.  The  flow  patterns  for  the  case  of  sinu.soidal  forcing  (symmetry 
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S  =  50%)  shown  in  Figure  2  are  similar  to  those  previously  observed  [4,5].  We 
note,  however,  that  the  shape  of  the  forcing  waveform  has  a  dramatic  influence  on 
the  flow  (see  Figure  3).  The  case  shown  in  Figure  3  illustrates  that  using  an  asym¬ 
metric  waveform  with  S  =  70%  results  in  a  larger  growth  rate  than  sinusoidal  forcing 
(Fig.  2,  4  Hz  case),  whereas  S  =  30%  produces  a  "forced"  layer  which  is  more  simi¬ 
lar  to  a  "natural"  layer. 

Figures  4,  5  demonstrate  the  effects  of  forcing  on  the  pdf  of  the  mixed-fluid 
concentration  field  P{^)  at  a  downstream  distance  of  about  x  =  1\  cm.  At  this  loca¬ 
tion,  the  Reynolds  number  of  the  natural  layer,  based  on  the  visual  thickness  and  the 
velocity  difference,  is  about  6,000  which  corresponds  to  a  turbulent  mixing  layer  dur¬ 
ing  the  mixing  transition.  The  most  striking  result  is  that  both  the  amount  and  the 
composition  of  mixed  fluid  are  modified  upon  external  forcing.  For  the  cases  we 
have  studied  thus  far,  the  predominant  concentration  shifts  to  larger  values  and  the 
total  amount  of  mixed  fluid  increases  (Fig.  4a).  Another  important  finding  is  that  the 
details  of  the  forcing  waveform  shape  have  a  significant  influence  on  the  mixing  field 
(Fig.  4b).  At  the  same  forcing  frequency,  changing  the  symmetry  of  the  waveform 
alters  both  the  amount  and  the  composition  of  the  mixed-fluid.  In  Figure  5,  we 
show  how  waveform  shaping  can  lead  to  either  major  (Fig.  5a)  or  minor  (Fig.  5b) 
changes  in  the  mixed-fluid  composition  field. 

The  changes  in  the  mixed-fluid  pdf  reponed  here  can  be  the  result  of 
modification  to  three  processes:  (1)  the  engulfment  of  the  free-stream  fluids  in  the 
first  rolled-up  vortex  (or  vortices),  (2)  the  overall,  2-D  entrainment  field  into  the 
layer  farther  downstream  beyond  the  first  roll-up  location,  and  (3)  development  of 
streamwise  vortices  and  3-D  small  scales.  We  are  in  the  process  of  isolating  these 
three  effects  by  making  concentration  measurements  in  the  early  stages  of  vortex  for¬ 
mation  and  by  studying  the  spanwise  structure  of  the  forced  flow  (both  visualization 
and  concentration  measurements). 
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2.  EFFECT  OF  CURVATURE  ON  MOLECULAR  MIXING 
IN  A  TWO-STREAM  SHEAR  LAYER 


In  ^  separate  smaller  effort,  we  are  studying  how  shear  layer  curvature  may 
affect  molecular  mixing.  The  extensive  body  of  knowledge  which  exists  on  turbulent 
flow  curvature  effects  concentrates  on  the  characterization  of  the  velocity  field  (see, 
for  exaiiiple.  Ref.  9)  whereas  much  less  is  known  about  the  mixing  field.  Besides 
the  inherent  interest  in  curvature  effects  on  fluid  mechanics  phenomena,  use  of  curva¬ 
ture  as  a  possible  tool  for  (molecular)  mixing  control  needs  to  be  investigated.  Cur¬ 
vature  could  change  the  mixing  field  by  modifying  the  overall  large  scale  (2-D) 
entrainment  into  the  flow  and  also  by  changes  in  the  3-D  structure  of  the  flow  due  to 
the  Taylor-Goertler  instability.  The  study  we  describe  below  outlines  the  portion  of 
our  work  on  the  effects  of  curvature  on  the  2-D  structure  of  the  shear  layer.  We  are 
planning  to  also  address  3-D  effects  which  may  be  the  most  important  mechnism 
through  which  curvature  modifies  the  mixing  field. 

We  investigated  the  effect  of  curvature  on  the  2-D  structure  of  a  shear  layer 
both  theoretically  and  experimentally.  The  inviscid,  linear,  parallel-flow  stability 
analysis  of  spatially  growing  disturbances  [10,11]  was  utilized  to  stuoy  the  instability 
characteristics  of  a  curved  shear  layer.  It  can  be  shown  that,  for  a  curved  layer,  the 
disturbance  eigenfunction,  <[),  satisfies  the  modified  Rayleigh  equation  given  by 


4)  -I-  — (J)  - 
R 


2  U  +  U/R  -  U/R^ 
+  - 


U  -  p/a 


4)  =  0 


where  (  )  corresponds  to  d  ldy‘-,  (  )  to  didy  and  U{y)  is  the  mean  velocity  profile. 
In  the  equation  above,  a  =  is  the  complex  non-dimensional  wave-number  and 

P  is  the  non-dimensional  frequency  which  is  taken  as  pure  real  for  the  present  spatial 
calculations.  The  nondimensional  radius  of  curvature  is  given  by  /?  =  p/5  where  p  is 
the  (dimensional)  radius  of  curvature  and  5  is  the  shear  layer  thickness.  A  positive 
value  of  R  corresponds  to  the  high-speed  stream  on  the  outside  whereas  a  negative 
value  implies  that  the  high-speed  stream  is  on  the  inside  of  the  curved  layer.  We 

use  a  mean  velocity  profile  which  incorporates  the  initial  wake  component  near  the 

splitter  plate.  The  form  of  the  velocity  profile  and  the  solution  method  for  obtaining 
the  values  of  unstable  a,  P,  and  the  corresponding  shape  of  the  eigenfunction  pCv) 

are  similar  to  those  described  in  Reference  1 1  and  will  not  be  repeated  here. 

The  effect  of  curvature  on  the  instability  characteristics  of  a  r  =  Ui/Ut  =  3 
shear  layer  (i.e.  X  =  (r-l)/(r-)-l)  =  0.5)  is  shown  in  Figures  6  and  7.  In  these 
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figures,  the  curvature  value  of  /?  =  100  is  large  enough  to  represent  the  limit  of  a 
straight  shear  layer  (i.e.  no  curvature).  Similar  to  previous  results  [11],  two  modes 
of  instability  are  found.  Mode  one  (the  shear  layer  mode)  leads  to  the  usual 
Kelvin-Helmholtz  shear  layer  rollup,  whereas  mode  2  (the  wake  mode)  rollup  pattern 
resembles  a  wake  flow.  We  note  that  the  effect  of  curvature  (positive  or  negative)  is 
minimal  even  at  a  curvature  as  small  as  10.  Similar  conclusion  is  reached  when  we 
examine  the  vorticity  eigenfunction  (i.e.  vorticity  mode  shape)  at  maximum 
amplification  (see  Figure  8). 

To  study  curvature  effects  experimentally,  the  shear  layer  facility  was  modified 
to  accommodate  a  curved  test  section.  Visualization  of  the  flow  patterns  for  natural 
and  forced  conditions  is  shown  in  Figure  9.  The  similarity  between  straight  and 
curved  shear  layers,  both  natural  and  forced,  supports  the  conclusions  reached  based 
on  the  sability  analysis. 

Our  results  to  date  indicate  that  curvature  may  not  play  any  major  role  in  modi¬ 
fying  the  initial  development  of  the  shear  layer  vortices  and  the  fluid  entrained  into 
them.  The  influence  of  curvature  on  streamwise  vorticity  is  known  to  be  very 
strong.  How  this  affects  the  molecular  mixing  field  remains  to  be  shown. 
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Figure  1.  Digital  image  time  sequence  of  the  mixing  field  for 
the  natural  shear  layer.  Flow  is  from  right  to  left  and  time 
increases  from  top  to  bottom.  Two  different  time  sequences  are 
shown  to  indicate  that  the  large  scale  organization  of  the  flow 
is  not  always  apparent. 
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Figure  2.  Digital  image  time  sequence  of  the  mixing  field  for 
sinusoidal  forcing.  The  shear  layer  natural  frequency  at  the 
splitter  plate  tip  is  about  Fo  »  27  Hz. 
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Figure  3.  Digital  image  time  sequence  of  the  mixing  field  for 
nonsinusoidal  forcing  at  F  =  4  Hz. 
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Figure  6.  Effect  of  curvature  on  the  amplification  rate  (-ai)  of  instabilit.' 
Modes  1  and  2  correspond  to  the  shear  layer  and  wake  modes,  respectively. 
Positive/negative  curvature  corresponds  to  the  high-speed  stream  on  the 
outside/inside  of  the  curved  layer.  The  straight  (infinite  R)  shear  layer 
limit  is  represented  by  the  curves  for  R  =  100. 
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